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SUMMARY 


Some  differences  between  the  mathematical  treatment  of  magnetic  cores  with  a  linear 
characteristic  and  those  with  rectangular  hysteresis  loops  are  reviewed.  Neeteson's 
approximation  to  the  flux  switching  characteristic  of  the  latter  is  applied  to  a  circuit 
with  the  primary  winding  energized  by  a  current  step,  and  the  secondary  winding  ter¬ 
minated  in  a  series  combination  of  inductance  and  resistance.  After  several  changes 
of  variable  the  resulting  differential  equation  reduces  to; 

^  +  +C2Sinu)  f  =  C3, 
dt 

where  u  is  a  function  of  the  flux  in  the  core ,  t  is  time,  and  the  C's  are  constants. 
By  assuming  that  sin  u  is  constant  over  successive  small  time  intervals,  this  equa¬ 
tion  is  solved  by  means  of  the  Laplace  transform.  An  auxiliary  differential  equation 
for  the  secondary  current  is  solved  by  utilizing  an  integrating  factor  and  numerical 
integration.  Curves  of  flux  and  current  versus  time  are  given  for  the  cases: 

Cf  =  1  ;  C2  =  0  ,  1  ,  2  ,  4  ;  C3  =  20  . 
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Section  1 
INTRODUCTION 

By  approximating  with  a  semicircle  the  flux  switching  characteristic  of  a  ferrite  core 
possessing  a  rectangular  hysteresis  loop,  Neeteson  (Ref.  5),  has  simplified  consider¬ 
ably  the  analysis  of  switching  circuits  containing  such  cores.  Generally, the  resulting 
differential  equations  are  nonlinear.  Neeteson  has  presented  solutions  for  circuits  in 
which  a  primary  winding  on  the  core  is  energized  by  a  step  current  generator,  and  a 
secondary  winding  is  terminated  by  resistance,  by  inductance,  and  by  a  parallel  com¬ 
bination  of  resistance  and  inductance.  One  is  led  quite  naturally  to  inquiring  into  the 
effect  of  terminating  the  secondary  winding  with  a  series  combination  of  resistance 
and  inductance.  It  will  be  shown  that,  after  certain  transformations,  the  resulting 
differential  equation  has  the  form 

$"<Ci  f  ‘  S 

dt 

The  flux  in  the  core  is  a  function  of  the  variable  u  . 

This  paper  is  concerned  principally  with  a  method  for  solving  numerically  the  fore¬ 
going  differential  equafion.  A  few  examples  have  been  considered.  The  solutions  are 
presented  in  the  form  of  graphs  for  the  variable  u  and  for  the  output  voltage  and  cur¬ 
rent,  which  are  derivable  from  u  by  simple  operations,  time  being  the  independent 
variable.  Before  delving  into  the  method  of  solution,  however,  it  is  desirable  to  pre¬ 
face  a  few  remarks  concerning  the  differences  in  analytical  treatments  between  cores 
with  linear  characteristics  and  cores  with  square  hysteresis  loops. 
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According  to  the  well-known  analogy  between  magnetic  and  electric  circuits,  the 
magneto-motive  force  in  ampere-tums  NI  (MKS  units),  the  reluctance  R  in  ampere- 
turns  per  weber  and  the  flux  (p  in  webers  are  related  by  the  equation  (Ref.  6). 

MMF  =  NI  =  R(p  .  (1) 

Neeteson  uses  the  symbol  (AT)  for  NI  ,  so  that,  following  his  notation,  Eq.  (1) 
becomes 


(AT)  =  R^.  (la) 

For  the  square-looped  core  being  switched,  the  pertinent  equation  relating  output  volt¬ 
age  (which  is  equal  to  d<p/dt  in  webers/sec)  and  the  ampere  turns  ( AT )  is 

^  =  (AT)  ■  V  ■  =  (AT)  •  r  •  ^1  -  ( 1  -  (2) 


where 

r  is  a  function  of  both  the  material  and  the  shape  of  the  core,  and  has  the 
dimensions  of  resistance  in  ohms, 

and;  i((p)  =  y/l  -  (1  -  <p/(P^)^  is  a  nondimensional  relation  approximating  thp 
flux  switching  characteristic  of  the  core  in  terms  of  the  ratio  of  the  flux  (j)  which 
has  been  switched  to  (p^  ,  the  latter  being  one-half  the  total  flux  capable  of  being 
switched. 

Two  figures,  taken  from  Neeteson,  will  aid  in  clarifying  the  ideas  involved.  Figure  1-1 
shows  an  idealized  hysteresis  loop  for  a  core,  and  serves  to  define  <p^  .  Figure  1-2 
shows  the  curve  of  a  function  derived  by  Lindsey  (Ref.  4)  which  Neeteson  approximated 
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by  the  foregoing  f{(p).  While  the  hysteresis  loop  of  Fig.  1-1  pertains  to  the  low- 
frequency  alternating-current  properties  of  the  core,  both  curves  of  Fig.  1-2  concern 
the  switching  properties  of  the  core  to  the  extent  described  by  Eq.  (2). 


<t>  =  24> 


Fig.  1-1  Idealised  Hysteresis  Curve  Illustrating 


the  Meaning  of  <p  and  .  (Taken 
from  Neeteson. ) 


Fig.  1-2  Approximation  of  Lindsey's  Curve  (drawn 
in  broken  line)  by  a  Semi- Circle.  (Taken 
from  Neeteson.) 
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The  important  point  is  that  in  Eq.  (1)  or  (la)  the  flux  is  proportional  to  the  ampere- 
turns,  whereas  in  Eq.  (2)  the  rate  of  change  of  flux  is  proportional  to  the  ampere- 
turns.  Since  d<^/dt  (in  webers/second)  is  equal  to  a  voltage,  Eq.  (2)  actually  states 
that  the  core  acts  like  a  nonlinear  resistor,  r  •  ]/ 1  -  (  \  .  in  contrast  to 

the  customary  consideration  of  a  winding  on  a  magnetic  material  as  an  inductor. 

The  magnetic  switching  resistance  r  ,  which  is  a  constant  multiplier  of  the  flux 
fimction,  depends  on  the  core  material  and  geometry.  Its  defining  equation  is: 

A 

r  =  ft  Y  ohms  ,  (core)  (3) 

where; 

p  is  the  magnetic  switching  resistivity  of  the  material  in  ohms  per  meter, 

f  is  the  mean  path  length  in  meters, 

A  is  the  cross-sectional  area  in  square  meters. 

This  equation  should  be  contrasted  with  the  well-known  formula  for  the  electrical 
resistance  of  a  wire  of  length  f  and  cross  sectional  area  A  : 

£ 

T  =  p  Y  ~  ohms  ,  (wire)  (4) 

where,  now,  p  is  the  electrical  resistivity  of  the  material  in  ohm-meters.  It  is 
immediately  evident  that  the  geometrical  terms  A  and  f  for  a  core  are  inverted  with 
respect  to  the  terms  for  a  wire;  consequently  the  magnetic  switching  resistivity  for  a 
core  has  the  dimension  ohms  per  meter.  It  is  interesting  to  see  why  this  inversion 
exists . 

The  core  resistance  is  defined  as  the  ratio  of  the  maximum  pulse  output  voltage 
(d(/>/dt)^  to  the  effective  input  ampere  turns  (AT): 
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Core  resistance  is  measured  experimentally  as  the  slope  of  the  straight  line  obtained 
by  plotting  the  maximum  pulse  output  voltage  (d(;p/dt)^  versus  the  input  ampere 
turns  1^^  ,  of  the  same  form  used  to  drive  the  circuit.  The  input  turns  1^^  are 
related  to  the  effective  input  ampere  turns  (AT)  by  the  relation 

(AT)  =  Njlj  -  (AT)q  ,  (6) 

where  (AT)^  is  the  threshold  value  below  which  no  output  voltage  due  to  irreversible 
switching  is  obtained.  See  Fig.  1-3  (taken  from  Neeteson).  Equations  (5)  and  (6)  may 
be  combined  to  give: 


(m)  =  -  r(AT) 

'  'm 


(7) 


Fig.  1-3  Graph  by  Means  of  Which  the  Value  of 

r  =  tan  o;  can  be  ascertained.  (Taken 
from  Neeteson.) 
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Upon  taking  the  derivative  of  ( d<^/dt  with  respect  to  ,  ( r  and  ( AT  are 

assumed  to  be  constant  with  respect  to  time)  we  obtain; 


i^) 

'  '  m 

d(NjIj)  • 


(8) 


This  expression  is  the  slope  of  the  foregoing  mentioned  straight  line.  The  quantities 
on  the  right-hand  side  of  Eq.  (8)  may  be  expressed  in  terms  of  the  more  fundamental 
quantities  B  and  H  (flux  density  and  field  intensity,  respectively): 

(f)  -(¥}  =*(f) 

'  'm  '  'm  '  'm 

and: 


With  these  substitutions  Eq.  (8)  becomes: 


,  r*/dB\  ./dB\ 

_  ■  ^  'nal  _  ^  m  A 

d(Hf)  dH  ■  f 


(9) 


The  derivative  d  |(dB/dt)^  j/dH 
quantity  has  been  denoted  by  p; 


is  a  function  of  the  core  material  only. 


P 


This 


(10) 


Equation  (9)  becomes: 


r 


ohms  . 


(core) 


The  foregoing  explanation,  beginning  with  Eq.  (5),  accounts  for  the  inversion  of 
geometrical  terms  referred  to  previously. 
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Equation  (2)  is  utilized  in  Section  2  to  derive  the  differential  equation  for  the  circuit 
of  concern.  This  is  followed  in  Section  3  by  the  mathematics  of  the  method  for  solving 
that  differential  equation.  In  Section  4  another  differential  equation  relating  the 
secondary  current  to  the  rate  of  change  of  flux  is  solved  by  means  of  an  integrating 
factor  and  numerical  integration.  Computed  curves  are  presented  and  discussed  in 
Section  5.  The  paper  is  concluded  with  a  short  summary  in  Section  6. 
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Section  2 


DERIVATION  OF  THE  DIFFERENTIAL  EQUATION 

.  ,  du  - 

-2  (Cl+  dt  =  ^3- 

dt 


The  circuit  to  be  analyzed  is  shown  schematically  in  Fig.  2-1.  Upon  taking  into  account 
the  secondary  ampere-turns  N2  I2  .  which  tend  to  oppose  the  driving  ampere-turns 
Nj  Ij^ ,  Eq.  (2)  of  Section  1  is  modified  to  become: 


d(p 

dt 


=  rl/l  -  (1  -  Kli  -  (AT)o  -  N2I2I 


(1) 


R 


2 


Fig.  2-1  Circuit  of  a  Core  Driven  by  a  Current  Generator  1^  , 
with  the  Secondary  Terminated  by  Inductance  and 
Resistance  in  Series. 
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In  order  to  find  the  secondary  current  we  make  use  of  Faraday's  law  that  the  voltage 
induced  in  the  secondary  winding  by  the  change  of  flux  in  the  core  is  Ng  d<^/dt  .  This 
is  equal  to  the  sum  of  the  voltage  drops  across  inductor  Lg  and  resistor  Rg : 


N 


d<p 
2  dt 


"2  dt 


(2) 


These  two  equations  suffice  to  define  the  circuit  behavior  completely.  The  driving 
force  is  assumed  to  be  a  known  fimction  of  time  as,  for  example,  a  step  function 

or  a  ramp.  The  maximum  flux  and  the  threshold  value  of  ampere  turns  (AT)q 

are  assumed  to  be  known,  as  well  as  the  values  of  Rg  ,  ,  and  N2  .  The 

quantities  (p  and  are  the  unknowns.  The  problem  consists  in  eliminating  one  of 
these,  Ig  I  for  example,  from  either  of  the  two  equations,  thereby  obtaining  an  equa¬ 
tion  in  terms  of  <p  only.  If  the  equation  in  <p  alone  can  be  solved,  the  answer  can  be 
substituted  into  Eq.  (2)  which  is  then  solved  for  . 


Proceeding  along  these  lines,  we  first  arrange  Eq.  (1)  in  the  form: 


'2  ■  "2 


"ill  -  (AT),,  - 


d(/> 


r  /l  -  (1  - 


(3) 


By  multiplying  both  sides  of  Eq.  (3)  by  R2 ,  and  by  multiplying  and  taking  the 
derivative  with  respect  to  time,  we  then  obtain: 


^2 

p  T  —  _ £ 

”•2^2  Ng 


-  (AT)o  - 


VTT 


(1  - 


2  dt 


(4) 


5.  h 

"2  dt  N„ 


^  d/ _ 1 


-  (d/w  ) 

'T'  Vjjj  / 


_  ^ 

2  dt 


) 


(5) 
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Substitution  of  Eqs.  (4)  and  (5)  into  Eq.  (2)  gives: 

dl. 


N  ^  .  3 

^2  dt  N2 


N 


"1  d 


1  ^^\ryll  -  (1  -  cp/<p^) 


d<p 
2  dt 


5 

N„ 


-  (AT),  -  ^  t 

rV  1  -  (1  -  <P/<P^) 


(6) 


After  similar  terms  are  grouped,  the  following  equation  results: 

h  d /  —  1.  —  dA ^  L  1 — =J\ 

N2rdt^^/7T77T;7,^)2dty  \^2  N^r 


2  1  dt 


d0 


N„ 


N^Il  -  (AT), 


+  L, 


'2  Ng  dt 


(7) 


Now,  let: 


by  means  of  which  we  obtain 


and: 


J£-  = 
“^m 


dw  _  „  dx 
dt 


Vl  -  (1  -  <p/v^)^  =  Vl  -  (1  -  x)^  . 
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With  these  substitutions  Eq.  (7)  becomes: 


Lo</’ 


’2'^m  d 


,  2'^m  _.i _ \ 

W  T  /  ;  2/ 

1  -  (1  -  x)7 


2'^m  ■  N„r  ^ 


dx 

dt 


N,  dl. 


-  ?r  K'l  -  'ATIol  *  S  TT  -dT 


(8) 


Next,  let: 


1  -  X  =  y 


whence: 


dx  „  dy 
'  dt  dt 


Substitution  of  these  into  Eq.  (8)  yields: 


jd  /  1  ,  _  L 

V  V 

^2  1  1  ^1 

=  (n,I,  -  (AT)„j  +  s;  dT 


R„  (p 
^  2 

*’m  N.r 


1  -  y^) 


dt 


(9) 


Lastly,  let: 


y  =  cos  u 
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whence: 


dt 


sin  u 


du 

dt 


1 


sin  u , 


and: 


A  /  1  =  A  /  sin  u  du^ 

dt  I  _  y2 


2 

d  u 
dt2 


With  this  new  change  of  variable  Eq.  (9)  is  transformed  into: 
L„ 


R. 


dr  ' 

'«2^m 

+  N»<p  sin  u 
2^m 

i 

N  dl 

Nili  - 

(AT)q 

"1^2  ^  dT 

du 

dt 


Upon  dividing  through  by  becomes 

•  ^  vr\  '  '  TYl 


d^ 

dt* 


"2'^m 


N,r  dL 


dt 


(10) 


(11) 
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For  the  particular  case  in  which  the  driving  force  is  a  step  function,  its  derivative 
with  respect  to  time  dlj^/dt  is  equal  to  0  for  t  5:  0  +  ,  whereupon  Eq.  (11)  reduces 
to: 


n2  ^ 

2  4 


(12) 


This  last  equation  has  the  form  indicated  in  the  introduction.  Before  presenting  its 
solution,  we  believe  it  worthwhile  to  review  the  transformations  leading  to  the  final 
form. 

Figures  2-2a,  2-2b,  and  2-2c  show  plots  of  the  function  t((p)  =  -  (1  - 

before  and  after  being  subjected  to  the  changes  of  variable:  x  =  fp/<P^  >  y  =  1  -  x  , 
The  arrows  on  the  curves  indicate  the  direction  taken  by  a  moving  point  as  (p 
increases.  In  Figs.  2-2a  and  2-2b  this  motion  is  clockwise,  while  in  Fig.  2-2c  this 
motion  is  counterclockwise.  The  contrast  of  directions  may  also  be  seen  by  tabulating 
corresponding  values  for  <p ,  x  and  y  as  follows: 

X  y 

0  1 

1  0 

2 


Now  consider  the  two  functions  y  =  cos  u ,  and  f(  u ) 
over  the  same  range  of  variables,  namely: 

y  =  cos  u  u  f(u 

1  0 

0  Tr/2 

- 1  tt 


lA  2 

=  yi  -  cos  u  =  sin  u, 


)  =  sin  u 
0 


<P 

0 

<P 

^m 

2<p 


m 
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Plots  of  these  functions  are  given  in  Figs.  2-3a  and  2-3b.  Again  the  arrows  on  the 

curves  indicate  the  direction  taken  by  a  moving  point  as  cp  increases .  Attention  is 

called  to  the  fact  that,  as  <p  increases  from  0  through  (p_  to  2<^„  ,  the  variable 

m  m 

u  increases  from  0  through  ■n/2  to  tt  . 

It  should  be  kept  in  mind  that  (p  is  the  dependent  variable,  whereas  time  is  the 
independent  variable.  The  values  <p  =  0  and  u  =  0  correspond  to  t  =  0,  and 
increasing  values  of  <p  and  u  are  associated  with  increasing  time.  When  t=  0, 
then  sin  u  =  0,  and  the  differential  Eq.  (12)  becomes: 


K'l- 


(13) 


At  the  time  when  cp  =  <p^  ,  sin  u  =  1 ,  and  the  differential  equation,  therefore,  has 
the  form: 


n: 


du  _ 
dt 


(14) 


At  the  time  when  cp  =  2(p^  ,  the  time  when  the  switching  process  ends,  we  again  have 
sin  u  =  0;  consequently,  the  differential  equation  once  more  reduces  to: 


^2  du 

S  dt 


Vm 


N^Il 


(15) 


It  is  not  difficult  to  find  an  electrical  analogue  for  Eq.  (12).  Assume  a  circuit 

^2^  r  ] 

consisting  of  a  battery  -  Pl^l  "  ®*'®^'Sizing  a  circuit  comprised  of  an 

^  m 

inductance  L2  in  series  with  a  fixed  resistor  and  a  variable  resistance.  Let 
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0 


(b) 


3  Plots  of  Functions  y  =  cos  u ,  and  f ( u )  =  -  cos^  u 
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the  variable  resistance  be  a  sinusoidal  function  of  the  charge  q  which  has  flowed 

2  2 

through  the  circuit  and  be  of  maximum  value  Ngr,  namely ,  Ngr  sin  q .  The  circuit 
is  shown  schematically  in  Fig.  2-4. 


Fig.  2-4  Electric  Circuit  Analogue  for  Eq.  (12) 


The  appropriate  electrical  differential  equation  is: 


n: 


r 


sin  q^ 


^  = 
dt 


(16) 


As  time  progresses  from  t  =  0  (at  which  time  q  =  0 ),  through  the  time  when 

q  =  7r/2  coulombs  to  the  time  at  which  q  =  tt  coulombs,  the  variable  resistance 
2 

ranges  from  0  to  Ng  r  and  back  to  0 . 

At  any  given  instant  of  time,  Eq.  (12)  may  be  considered  as  a  second  order  differential 
equation  with  constant  coefficients.  At  a  slightly  later  instant  of  time  Eq.  (12)  may 
again  be  considered  to  have  constant  coefficients,  but  with  the  coefficient  of  the  first 
derivative  a  little  altered  from  its  previous  value.  This  suggests  that  a  possible  method 
of  solution  is  to  solve  the  differential  equation  for  each  particular  instant  of  time  by  a 
standard  method  as,  for  example,  by  the  Laplace  transform,  using  as  the  initial  con¬ 
ditions  for  each  solution  the  value  of  the  dependent  variable  and  its  first  derivative 
obtained  from  the  solution  for  a  slightly  prior  time.  This  solution  was  completed  and 
will  be  outlined  in  Section  3,  For  that  solution  we  need  first  to  settle  ig)on  a  value  for 
du/dt  at  t  =  0 ,  which  will  occupy  us  for  the  rest  of  Section  2.  (The  value  u  =  0  at 
t  =  0  is  the  second  initial  condition  needed,  but  this  is  already  known. ) 
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Lacking  any  definite  knowledge  about  du/dt  at  t  =  0 ,  let  us  assume  that  it  is  equal 
to  zero.  We  shall  now  show  that  this  assumption  implies  that  in  the  vicinity  of  t  =  0 
the  flux  (p  must  vary  as  a  power  of  t  greater  than  2  .  Since  in  view  of  experimental 
work  this  is  reasonable,  the  assumption  will  be  used  in  the  later  treatment. 

From  the  relation; 


u  =  arc  cos  y 


we  get,  by  differentiation: 


du 

dt 


1  ^  _ _ 1  ^ 

V^i  -  y2  V^i  -  y  +  y 


(17) 


At  t  =  0 ,  y  =  1  ,  and  the  denominator  of  Eq.  (17)  is  zero.  In  order  for  du/dt  to 


remain  finite  at  this  point,  we  must  have  ^  0  a  greater  power  of  (1-y)  than  ^  , 


say  as  (1-y)  where,  for  reasons  which  will  appear  presently,  the  range  of 


is  restricted  to  0  <  e  <  9  •  From  the  original  transformations: 


1  -  y  = 


X  = 

<p 

^  TV-k 


the  assumption  that; 


(1  -  y) 


1  + 


e 
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requires  that: 


Therefore  let 


d<p 

dt 


—  0 


as 


J2_ 

(P 


m 


d<p 

dt 


=  a  constant) 


or  alternatively: 


Kj  dt 


An  indefinite  integral  yields: 

1 

2  ~  ^ 

<p  =  Kjt  +  Kg  ( Kg  =  a  constant) 


(18) 


(19) 


(20) 


whence: 


<P  =  (K,t  + 


Kg) 


=  (Kjt  + 


^2) 


2/(l-2e) 


(21) 


When  t  =  0,  <p  =  0 ;  therefore.  Kg  =  0.  Consequently,  we  get: 

2 

^  =  K^t^"^^  .  (22) 


Since  e  is  positive,  <p  varies  as  a  power  of  t  greater  than  2  in  the  vicinity  of 
t  =  0. 
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K  e  =  1/2  ,  we  must  go  back  to  Eq.  (19)  and  repeat  the  integration.  We  then  obtain: 

^  =  K,  dt  (23) 

(pi  ' 

Imp  =  Kjt  +  an  Kg  (24) 


whence: 

^  =  Kg  exp  (K^t)  .  (25) 

The  value  t  =  0  yields  <p  =  K^,  which  is  inconsistent  with  the  previously  assigned 
condition  cp  =  0 .  Therefore,  the  value  e  =  1/2  is  not  allowed. 

For  values  of  e  >  1/2,  the  exponent  of  (Kj^  t  +  Kg)  in  Eq.  (21)  is  negative  and  the 
following  results: 


(p  = 


(K^t  +  Kg) 


2e-l 


The  value  t  ==  0  yields 


<P 


K, 


2€-1 


which  again  is  inconsistent  with  the  assigned  condition  (p  =  0  .  Therefore,  the  values 
e  >  1/2  also  are  not  allowed.  Consequently,  the  range  of  e  compatible  with  the  con¬ 
ditions  cp  =  0 ,  d<p/dt  =0  at  t  =  0  (or  u  =  0 ,  du/dt  =0  at  t  =  0 )  is: 


0  >  €  >  2  , 
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as  mentioned  above.  The  value  of  e  will  not  be  determined  precisely,  although  it 
could  be  determined  from  a  particular  solution.  The  discussion  on  the  range  of  e 
merely  serves  to  give  us  an  idea  of  the  shape  of  the  flux  curve  near  the  time  origin, 
that  is,  it  is  at  least  parabolic. 
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Section  3 

ITERATIVE  PROCEDURE  FOR  SOLVING  THE  DIFFERENTIAL  EQUATION 


For  reference  we  repeat  Eq.  (12)  of  Section  2: 


d 

dt 


^  u  /^2  ^2  r  \  du  ^2  r  y  /  a  t  \  1 


(i> 


For  brevity  we  shall  make  the  further  substitutions : 


A  2  ,  „2  r  , 

^  =  L^'"N2L^sinu(t.) 


B 


In  these  substitutions,  B  is  a  constant.  The  subscript  i  of  Aj  and  8lnu(tj)  is  used 
to  indicate  that  these  quantities,  which  are  variable,  wttl  be  evaluated  at  the  discrete 
times  tj(i  =  0,  1,  2,  . . . ),  and  considered  constant  over  the  interval  ^ -  t^  .  In 
addition  to  u(t.)  we  must  evaluate  its  first  derivative  with  respect  to  time  u'(tj). 

With  this  new  notation  Eq.  (1)  becomes: 


+  A 


du 

idt 


=  B 


(2) 


subject  to  the  initial  conditions 


u(tQ)  =  0 
u'(tQ)  =  0  . 
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i 


Since  we  are  considering  Aj  as  a  constant  for  the  present,  we  may  apply  the  Laplace 
transform  to  Eq.  (2).  To  do  this  we  need  the  following  relations  (Ref.  1): 


/j  d  u 

*s  2 

=  s^U(s)  -  su(t.)  -  U'(tj) 

(3a) 

* 

dt 

i 

/}  du 
<  dt 

=  sU(s)  -  u(t.) 

(3b) 

^  B 
s  ’ 

(3c) 

where  U(s)  is  the  Laplace  transform  of  u(t)  .  The  use  of  u(t.)  and  u'(tj)  in 
Eq.  (3)  rather  than  0(1^)  or  u'(tQ)  [or  u(0+)  and  u'(0+)  in  more  customary 
notation]  serves  to  emphasize  the  fact  that  time  is  measured  from  tj  (i=0,  1,  2,  ...) 
rather  than  only  from  t^  .  Substitution  of  Eq,  (3)  into  Eq.  (2)  yields : 

s2u(s)  -  su(t.)  -  u'(t.)  +  Aj  [sU(8)  -  u(t.)]  =  (4) 


Rearranging  the  terms  of  Eq.  (4)  we  obtain : 


whence: 


|s^  +  A.s|u(s)  =  ®  +  u'(t.)  +  A.u(t.)  +  su(t.) 


U(s)  = 


B 


u'(t.)  +  A.u(t  )  u(t  ) 
+ - = - = - =-  + 


s  (s  +  Aj) 


s  ( s  +  Aj )  s  +  Aj 


(5) 
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From  a  table  of  Lafdace  transforms  (Ref.  2)  we  find  the  following: 


0-1  1 


exp  (-Ajt)  +Ajt  -  1 


s  (s  +  A^) 


(7a) 


1  _  1  -  exp(-A.t) 

<  s(s  +  A.)  A. 


(7b) 


1 


(7c) 


With  the  aid  of  these,  the  inverse  transform  of  Eq.  (6)  is: 


u(t)  =  |exp  (-Ajt)  +  A.t  -  1  + 


u'(tj)  +  A.u(t.) 


[l  -  exp  (  -A.t)| 


+  u  ( tj )  exp  ( -  Aj  t ) 


(8) 


Multiplying  out  the  second  group  of  terms  on  the  right  and  combining  it  with  the  third  we 
obtain: 


^(tj)f  iRf  1 

u(t)  =  u(t.)  +  -J—  jl  -  exp  (-A.t) j  +  -g  Jexp  (-Ajt)  +  A.t  -  l|  (9) 
1  A . 


Differentiation  of  Eq.  (9)  yields : 


u'(t)  =  u'(t.)  exp  (-Ajt)  +  ^  [l  -  exp  (-Ajt)j 


(10) 
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The  forms  of  these  equations  are  interesting.  Let  us  first  consider  Eq.  (9).  This 
equation  states  that  the  value  of  u  at  some  time  t  measured  from  time  t^  ,  namely 
u(t)  ,  depends  on  the  value  of  u  at  the  prior  time  u(t^)  ;  on  its  derivative  at  that 
time  u'(t. )  multiplied  by  [1  -  exp  (-A.t)]/A.  ;  and  on  the  driving  force  B  multiplied 

1  Q  1  1 

by  [exp  (-Ajt)  +  A.t  -  1  ]/A^  •  Equation  (10)  says  that  the  derivative  of  u  at  some 
time  t  measured  from  time  tj  ,  namely  u'(t)  ,  depends  on  the  value  of  u'  at  the 
prior  time  u'(t.)  ,  multiplied  by  exp(-Ajt)  ,  and  on  the  driving  force  B  multiplied 
by  [  1  -  exp  (  -  A.  t )  1  /A .  .  These  equations  are  assumed  to  be  valid  only  over 
a  small  time  interval  At  ~  +  i  •  These  equations,  therefore,  should  more  prop¬ 

erly  be  written  as 


u'(t,) 


“'•i  +  i*  '  "<'i>  ♦  -IT  I'  -  ('Viti  ■  ‘i 


+  -l,)|  (11) 


and; 


u'(t,,il 


u'(t,)exp(-Aj(t,^j  -  t.))  [l  -  exp|-A,(t,^j  -  t,)])  (12) 


From  the  foregoing  equations ,  it  is  apparent  that  if  we  know  the  values  of  u ,  u'  and  A 
at  time  t.  ,  namely  u(t.)  ,  u'(t.)  ,  and  A.  ,  we  can  then  use  Eqs.  (11),  (12)  and  the 
following  relation: 


to  calculate  u(tj^j)  ,  u'(t.  ^j^)  and  .  These  last  equations  may  then  be  utilized 

to  calculate  u(t.  ^2)  »  '^i  +  2  ’  Proceeding  in  this  way,  we  are  able 

to  calculate  the  solution  of  Eq.  (1),  starting  with  u(tQ)  =  0  ,  and  continuing  until 
u(t.)  =  TT. 
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This  procedure  lends  itself  readily  to  machine  or  hand  computation.  It  will  be  noticed 
that  the  expression  multipling  u'(t. )  in  Eq.  (11)  is  identical  with  that  multiplying  B 
in  Eq.  (12).  Consequently,  only  three  exponential  expressions  need  be  evaluated  at 
each  step. 


If  in  Eqs.  (11)  and  (12)  we  denote  the  time  interval  by  At  “  choose 

this  interval  small  enough  so  that  A^At  «  1  ,  then  it  is  convenient  to  expand  the  expo¬ 
nential  expressions  in  infinite  series  of  which  only  the  first  few  terms  will  be  needed. 
We  then  obtain: 


r  b  -  e^p(-Vt))=  ^ 


1  - 


(A  At)^  (A  At)^ 

1  - ^  + 


A. 


=  At 


(AAt)^  (AAt)^ 
■^i^*  2!  '''  3l 


A. At  (A. At) 
1 - ^  +  * 


2!  31 


(14a) 


^  |exp(-A.At)  +  A. At  -  l] 


a: 

1 


A^ 

1 


(A  At)^  (A  At)^  (A  At)*^ 


(A.At)^  (A  At)^  (A  At)^ 
1  1 


21 


31 


41 


+  AjAt  -  1 


=  (At)' 


j  A.  At  (AjAt) 

^  W  41 


(14b) 
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exp  (-AjAt)  =  1  A.  At  +  — ^ 


(AjAt)^  (AjAt)^ 


31 


(14c) 


With  these  substitutions,  Eqs.  (11)  and  (12)  become: 


«(ti  +  i)  =  u(t.)  +  u'(t.)At 


A  At  (A  At)^  (A.At)^ 

1  -  — —  +  ^  ^  . 

2!  ^  3!  “  41 


+  B(At)' 


1  A.At  (A.At)2  (A.At)^ 


2!  3! 


41 


51 


(15) 


and: 


“'(ti  +  l)  =  u'(t^) 


(A  At)^  (A.At)^ 

^  ^  - Vr-  - 


+  BAt 


A  At  (A  At)2  (A.At)^ 

1 _ i 1  t 

21+31  ^1  + 


(16) 


As  an  example  of  the  use  of  these  formulas,  consider  the  following  equation: 


+  ( 1  +  sin  u)  ^  =  20 
dt 


(17) 


Here  we  have  set: 


2 

T—  =  1  second 
^2 
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N?  ^  =  1  second  ^ 


B  =  7^  —  [n,  I,  -  (AT).]  =  20  seconds 

4  I  1  1  OJ 


Then: 


A.  =  l+  sinu(tj)  (i=0,l,2,...)  . 

As  u  varies  from  0  through  7r/2  to  tt  ,  the  quantity  A.  varies  from  1  through  2 
and  back  to  1  again.  Let  us  choose  a  time  interval  At  =  0.05  seconds  ,  which  is 
small  compared  to  the  time  constants  ^>2/^2  ^2^^2  ^  which  is  equal  to 

1  second.  We  also  note  that  A.  At  varies  between  0.  05  and  0. 1  which  are  small  com¬ 
pared  to  1,  so  that  only  a  few  terms  of  the  expansions  in  Eqs.  (15)  and  (16)  are  required. 
Utilizing  the  conditions 


u(to)  =  0 


u’(tQ)  =  0 


we  have  that : 


A^At  =  0.05 

and  that: 

u(tj)  =  u(0.05)  =  0  +  0  +  20(0.05)^  (0.491  770)  =  0.024  589 
u'(t^)  =  u'(0.05)  =  0  +  20(0. 05)(0. 975  412)  =  0.975  412 
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sinu(tj)  =  sin  (0.024589)  =  0.024  587 
Aj  =  1  0.024  587  =  1.024  587 

For  the  next  set  of  computations  we  then  get; 

AjAt  =  0.051  229 

uftg)  =  u(0. 10)  =  0.024  589  +  0.975  412 ( 0.  05)( 0.  974  817) 

-H  20(0.05)^(0.49  1  570)  =  0.09  6  710 

u'(t2)  =  u'(0.  10)  =  0.975  412(0.950  061)  -i-  20  ( 0.  05 )( 0.  974  817)  =  1.901  518 
and  so  forth. 

From  the  solutions  for  u  and  u'  we  are  then  able  to  determine  the  flux  and  its  rate 
of  change  by  going  back  through  the  transformations  : 

1  -  y  =  X  =  ,  y  =  cos  u  . 

m 

There  results: 

=  1  -  cos  u  (18) 

and: 

=  sinu^  =  sinu[u'(t)] 

in  which  w  is  assumed  to  be  known, 
m 
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Section  4 

SOLUTION  OF  THE  DIFFERENTIAL  EQUATIONS  FOR  THE  SECONDARY  CURRENT 


In  order  to  determine  the  secondary  current  we  have  to  solve  Eq.  (2)  of  Section  2  which 
is  rewritten  in  the  form ; 


R«  N„  , 

+  -2  T  =  ^  ^ 

^2  2  L^dt  - 


(1) 


Here  ^  is  known  from  Section  3,  but  not  in  a  neat  analytical  form.  Equation  (1)  is  a 
particular  case  of  a  more  general  equation  whose  solution  is  readily  available  (Ref, 7): 


dig 

^  +  P(t)l2  =  Q(t)  (2) 


Upon  multiplying  both  sides  of  this  equation  by  the  integrating  factor  exp  [  /  P(t)dt]  , 
the  left-hand  side  becomes  an  exact  differential  and  results  in  the  following: 

^[igcxp  [/P(t)dtl|  =  Q(t)  exp  [/P(t)dt]  .  (3) 

Integration  of  Eq.  (3)  then  gives ; 

Ig  exp  [/P(t)dt]  =  jQ(t)  exp  [/P(t)dt]dt  +  C  (4) 

where  C  is  a  constant. 
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Applying  this  method  to  Eq.  (1)  with: 


P(t) 


Q(t) 


^2 


we  obtain  the  counterpart  to  Eq.  (3) : 


dt 


(5) 


In  order  to  avoid  the  integration  constant  C  which  appeared  in  Eq.  (4) ,  we  shall  inte¬ 
grate  Eq.  (5)  between  the  limits  0  and  t  ,  rather  than  use  an  indefinite  integral.  We 
then  obtain: 


from  which: 


(6) 


(7) 


where  the  notation  l2(t)  is  used  to  emphasize  that  is  a  function  of  time ,  and  l2(0) 
is  the  current  I2  at  t  =  0  .  Solving  for  l2(t)  we  then  find: 


yt)  =  yo)  exp 


(8) 
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Since  in  our  case  12(0)  =  0  ,  Eq.  (8)  reduces  to; 


yt)  = 


(9) 


It  will  be  more  convenient  to  obtain  the  solution  in  terms  of  ^  ,  which  is 

□t  Qt 

easily  done  as  follows ; 


(P 


m 


dt 


(10) 


This  foregoing  equation  may  be  normalized  by  dividing  through  by  N„(/)  /L„  to  give: 

M  in '  dt 

(-1^^)  /  aF®^(Lf 


l2(t) 


No*;*’  /Lo 
2^m'  2 


=  exp 


The  integral  on  the  right  in  Eq.  (11)  was  evaluated  by  using  the  trapezoidal  rule  (Ref. 
7)  which  is  simple  to  apply  and  well  suited  to  provide  a  "running"  value  of  the  integral 
Multiplying  this  running  value  of  the  integral  by  exp  [-(R2/L2)  t]  gives  the  normal¬ 
ized  current.  By  way  of  illustration,  for  the  differential  equation  (17)  used  as  an  ex¬ 
ample  in  Section  3,  the  right-hand  side  of  Eq.  (11)  becomes: 


-t 

e 


J  dt 
0 


t 


e 


dt  . 
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dx 


The  first  few  values  of  arrived  at  by  the  iterative  procedure  of  Section  3 ,  and  of 


ample,  that  time  is  approximately  0.  65  seconds,  at  which  time  the  switching  process 

’m 


is  complete,  the  flux  cp  has  reached  the  value  ,  and  thereafter  ^  0  . 


dx  .t 
dt  ®  ’ 

dx  t 

J  dt  ’ 

0 

,  -t  1  dx  t,. 
and  e  y  ^  e  dt 

are  as  follows : 

t 

dx 

dt 

dx  t 
dt  ® 

f^dx  tj. 

-t  dx  t  .. 

e  /  e  dt 

i 

k 

0 

0 

0 

0 

0 

I 

0.05 

0.023  982 

0.025  212 

0.000  630 

0.000  599 

0.10 

0.183  609 

0.202  919 

0.006  334 

0.005  731 

0.15 

0.587  380 

0. 682  438 

0.028  468 

0.024  503 

The  integration  is  carried  out  from  t  =  0  to  the  time  at  which  u 

=  TT  .  For  this  ex- 

The  current  in  the  secondary  circuit  has  some  finite  value,  and  must  decay  to  zero 
through  the  inductance  ,  the  resistance  Rg  ,  and  the  secondary  turns  Ng  linking 
the  core.  Since  ^  is  zero,  and  will  be  assumed  as  remaining  zero  throughout  the 
current  decay  time,  there  is  no  voltage  induced  in  the  secondary.  The  secondary 
winding,  consequently,  may  be  assumed  as  a  short  circuit  during  this  time.  The  cur¬ 
rent  decays,  therefore,  from  its  previous  final  value  12^  in  accordance  with  the  differ¬ 
ential  equation: 


^2  dt  ^2^2 


(12) 
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The  solution  is  known  to  be  (Ref.  3): 


l2(t)  = 


(13) 


where  time  t  is  now  measured  from  the  time  at  which  the  switching  process  reached 
completion.  In  order  to  keep  the  quantities  involved  consistent  with  the  prior  notation, 
this  equation  may  be  normalized  also  and  written  as: 


l2(t) 


2 


(14) 


For  the  example  we  are  pursuing,  the  final  normalized  current  was  1. 595  at  t  =  0.  65 
seconds.  Thereafter  the  current  was  assumed  to  decay  as: 


l2(t) 


=  1.595  e 


-t 
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Section  5 

COMPUTED  CURVES 


The  following  differential  equations  were  solved  by  the  methods  outlined  previously. 


(1) 


+sinu)f^  20 
dt 


(2) 


^  +  (1  +  2  sin  u)  ^  =  20 


dt 


(3) 


^  +  (1  +  4  sinu)  ^  =  20 
dt^ 


(4) 


Curves  were  plotted  for  the  quantities  u(t)  ,  u'(t)  ,  x(t)  ,  x'(t)  and  l2(t).  with 
I2  ( t )  presented  on  two  different  time  scales  in  order  to  show  the  large  trailing  edge 
of  the  current.  These  are  displayed  for  Eq.  (1)  in  Figs.  5-la,  5-lb,  5-lc;  for  Eq.  (2) 
in  Figs.  5-2a,  5-2b,  5-2c;  for  Eq.  (3)  in  Figs.  5-3a,  5-3b,  5-3c;  and  for  Eq.  (4)  in 
Figs.  5-4a,  5-4b,  5-4c.  The  ordinate  and  abscissa  corresponding  to  the  time  when 
the  variable  u  =  tt  ,  which  ends  the  switching  process,  are  indicated  on  Figs.  5-la, 
5-2a,  5-3a,  and  5 -4a.  As  the  coefficient  of  sin  u  in  the  differential  equations  varies 
from  0  through  4  ,  the  switching  time  increased  from  0.  615  through  0.  655  and  0.705 
to  0.830  seconds. 

The  reader  is  probably  most  familiar  with  the  curve  of  u'  for  Eq.  (1) ,  which  corres¬ 
ponds  to  the  rise  of  electric  current  in  a  series  L-R  circuit  with  fixed  circuit  elements. 
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t  -  seconds 


t  -  seconds 


Fig.  5-lc 


Plot  of  + 

dr 


du 

dt 


20:  Normalized  Current  vs.  t 
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t  -  seconds 


,2  , 

Fig.  5-2a  Plot  of  +  ( i  +  sin  u )  ^  =  20 ;  u  and 

dt 
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0  .2  .4  .6  .8  1.0  1.2  1.4  1.6  1.8  2.0  2.2  2.4  2.6  2.8  3.0  3.2  3.4  3.6 

t  -  seconds  _ 


Fig.  5-2c 


,2 

Plot  of  — ^  +  (1  +  sin  u)  ^  =  20: 
dt 


Normalized  Current  vs.  t 
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i 


0  .1  .2  .3  .4  .5  .6  .7  .8  .9 

t  -  se  conds 


I 

I 


1 
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t  -  seconds 


Fig. 


i-4b  Plot  of  +  ( 1  +  4  sin  u )  ^ 

dr 


20  :  X  and  x'  vs.  t 
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energized  by  a  battery.  In  view  of  the  analogy  cited  In  Section  2  (see  Fig.  2-4)  the 
pertinent  equation  for  the  circuit  charge  Is: 

^.^.20  ,5, 

dt 

If  we  change  the  dependent  variable  from  charge  to  current  by  means  of  the  following 
substitution: 


I 


=  ^ 
dt 


(6) 


this  becomes  : 


20 


(7) 


The  solution  Is : 


I  =  20(1  -  e"S  (8) 

which  represents  an  exponential  rise  of  current  to  a  final  value  of  20  amperes  with  a 
time  constant  of  1.  The  variable  u'  has  the  following  identical  form: 

u'  =  20(1  -  e"S  (9) 

of  which  only  8/10  of  a  time  constant  appears  in  Fig.  5-la.  Since  the  curve  for  u'  is 
almost  linear  over  this  portion,  the  curve  for  u  is  almost  parabolic.  As  the  coeffi- 
dent  Ng  “  increases  (corresponding  to  increasing  Ng  or  r  or  both, 

L„  remaining  fixed) ,  a  pronounced  wiggle  appears  In  the  curve  for  u'  .  This  finally 

A 

tends  to  make  u'  vary  about  a  constant  value  over  a  portion  of  its  range  (as  u'  = 
4.375  ±  .  125  over  the  range  0.  285  s  t  ^  .725  in  Fig.  5-4a)  which  causes  the  variation 
of  u  to  change  from  parabolic  to  almost  linear. 
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The  output  voltage ,  which  is  proportional  to  x'  ,  has  the  highest  peak  for  the  coeffici¬ 
ent  of  sin  u  equal  to  zero,  then  decreasing  from  7. 1  through  6.3  and  5.  6  to  4.45  for 
a  coefficient  of  4.  Peak  values  of  corresponding  normalized  output  currents  decrease 
from  1.  64  through  1. 59  and  1.  56  to  1. 48. 
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Section  6 
CONCLUSIONS 


Starting  from  Neeteson's  approximation  for  the  flux  switching  characteristic  of  a  mag¬ 
netic  core,  we  have  derived  the  differential  equation  describing  the  case  when  such  a 
core  is  terminated  in  series  inductance  and  resistance.  By  means  of  successive  trans¬ 
formations  ,  the  equation  was  reduced  to  a  form  which  could  be  solved  by  an  iterative 
procedure.  A  family  of  solutions  was  computed  and  plotted. 

To  date  no  attempt  has  been  made  to  check  these  curves  against  experimental  results 
although  verification  is  planned  for  the  future.  Since  Neeteson's  approximation  is  ad¬ 
mittedly  a  compromise  between  exact  fit  to  the  flux  characteristic  and  mathematical 
simplicity,  only  a  fair  agreement  is  to  be  expected. 

It  is  believed  that  the  differential  equation  and  solutions  are  applicable  to  flux  steering 
through  multipath  magnetic  structures. 

The  mathematical  methods  employed  are  themselves  of  some  interest. 
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